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Abstract. When the Lake Louise Winter Institute started twenty-five years ago, many 

I I properties of quantum chromodynamics (QCD) were believed to be true, but had not been 

,S^ demonstrated to be true. This talk surveys a variety of results that have been established with 

Mh lattice gauge theory, directly from the QCD Lagrangian, shedding light on the origin of (your) 

mass and its interplay with dynamical symmetry breaking, as well as some further intriguing 

features of the natural world. 



^ 1. Solving QCD 

■^ Quantum chromodynamics (QCD) is the modern theory of the strong nuclear force. It is part of 

'^ the Standard Model of elementary particles, yet also has profound influence on nuclear physics 

and on astrophysics. It is also rich and fascinating. In this talk, I aim to cover some results that 
are interesting in their own right, influential in a wider arena, quantitatively impressive, and/or 

^^ qualtitatively noteworthy. 

T— I The Lagrangian of QCD has "1 + n/ + 1" free parameters: 

"^ rif - 

X ^QCD = ^ tr[F^,F^n - Y^ i,0 + ms)^s + ^.e^"''' tr[F^,F,<,], (1) 

where F^'^ is the gluon's field strength, p = 7^(9^ + A'^), and -0/ denotes the quark field of 
flavor /. The first parameter is the gauge coupling g^ , the next nj are the quark masses nif, 
and the last, 6, multiplies an interaction that violates CP symmetry. There are six quarks (that 
we know about), but at energies below the top, bottom, and charm thresholds, it is convenient 
and customary to absorb the short-distance effects of these quarks into a shift of g^ and then 
take QCD with re/ = 5, 4, or 3. In addition to these shifts, the coupling (7^ diminishes gradually 
with increasing energy, stemming from virtual processes of gluons and the Uf active quarks; 
this is called "asymptotic freedom" [ll [2] . More generally, one could imagine a matrix in the 
mass term, V'aiT-abV'fe) with eigenvalues mf in Eq. n\\. In this context, the coupling multiplying 
e^^P" tT:[Ffj_yFpfj\ is altered: 6 = — argdetm. In the Standard Model, 9 is considered purely 
chromodynamic, whereas ruab arises from Yukawa couplings between quarks and the weak- 
isodoublet Higgs boson. Only the difference 9 is observable. 

Before saying that a mathematical theory describes or explains the natural world, one must 
fix the free parameters with the corresponding number of measurements, in this case 1 + rej + 1. 
Because the color of quarks and gluons is confined, the free parameters of QCD must be 
connected to properties of QCD's eigenstates, which are the bound states called hadrons. At 



high energies, a sum over many hadronic states can be related to a sum over many quark- 
antiquark-gluon states ( "quark-hadron duahty" [5]). At lower energies, this is not possible. In 
this nonperturbative regime, it is preferable to relate all the parameters to quantities like hadron 
masses, whose experimental interpretation is clear. Then the obstacle is (merely) to compute 
the relationship between the QCD Lagrangian and hadronic properties. 

A long-promising and now-successful approach to such computations is to formulate QCD 
as a lattice gauge theory [4j. Then the ultraviolet cutoff (needed in any calculation) is built in 
from the outset, and the correlation functions of QCD are mathematically well-defined. The 
parameters are fixed as follows. The electric-dipole moment of the neutron is unobservably 
small, leading to a bound 9 < 10~^^. Such delicate cancellation of 6 and argdet ?7i is a mystery, 
known as the strong CP problem |5j, but for QCD calculations it simply means we can set 
^ = with no important consequences. The rest are tuned to reproduce 1 + ?^/ specific hadronic 
properties. Even the dimensionless gauge coupling is related to a dimensionful, measurable 
quantity, because the gauge coupling runs. Thus, to lend a physical interpretation to g'^ one 
has to compute the energy at which the coupling reaches a fiducial value, like g^ = 1. But the 
calculations don't know a priori about kg or GeV/c^; instead the energy at which g^ = 1 can 
be computed only relative to some other standard mass, such as the mass of the nucleonjj 

It is worth solving QCD quantitatively for many practical reasons in particle physics, nuclear 
physics, and astrophysics. In such applications, the enterprise is unabashedly pragmatic, with a 
focus on a clear understanding of uncertainties remaining in the calculation. This is important 
but a secondary focus of this talk. 

The primary focus is as follows. The conception of QCD is rightly hailed as a triumph of 
reductionism, melding the quark model, the idea of color, and the parton model into a dynamical 
quantum field theory. The application of QCD is, however, rich in emergent phenomena. 
Symmetries emerge in idealized limits: C, P, and T are exact when ^ = 0; chiral symmetries 
emerge when two or more quark masses vanish [6]; and heavy-quark symmetries are revealed 
as a quark mass goes to infinity [71 Ej^. More remarkable still are the dynamical phenomena 
that emerge, starting with Aqcc the "typical scale of QCD," which is not an input. Much of 
what is "known" about QCD in this essentially nonperturbative arena has been, for a long time, 
based on belief: Evidence from high-energy scattering led to the opinion that QCD explains all 
of the strong interactions. This opinion led to the belief that QCD exhibits certain properties, 
because otherwise it would not be consistent with observations. These emergent phenomena — 
such as chiral symmetry breaking, the generation of large hadron masses despite very small 
quark masses, and the thermodynamic phase structure — are the most profound phenomena of 
gauge theories. The aim of this talk is to survey how lattice QCD has filled many of these gaps, 
replacing belief with knowledge. 

The rest of this talk is organized as follows. Section [2] give a short review of lattice-QCD 
methodology and jargon. Hadron masses and their connection to chiral symmetry are discussed 
in Sees, p^ and Hi An output of these calculations are the (in some cases remarkably small) quark 
masses and the gauge coupling; these results are discussed in Sec. [5| along with a few "tense" 
results on flavor physics. The phase structure of QCD is discussed in Sec. [6j Section [7] offers 
some perspective. The developments reported here took place during quarter-century history 
of the Lake Louise Winter Institute. Many of the leading players have been Canadian, and 
wherever I know of a Canadian connection I've marked the contribution with a maple leaf W. 

2. Lattice Gauge Theory 

Lattice gauge theory \4\ was invented in an attempt to understand asymptotic freedom without 
introducing gauge- fixing and ghosts [9]. The key innovation of Ref. [4J is to formulate non- 

^ In principle, the nucleon mass is a good example, but it depends sensitively enough on quark masses that, in 
practice, one uses strange baryon masses, quarkonium mass splittings, or other quantities that are not as sensitive. 



Abelian gauge invariance on a spacetime lattice. Then the functional integrals defining QCD 
correlation functions are well-defined: 

(•) = ^ y VAVi^Vi; [.] exp {-S) , (2) 

because the measures "DU, "Dtp, Dip are products of a countable number of normal differentials. 
Here S = f d^xC is the action, • is just about anything, and Z ensures (1) = 1. This formulation 
is formally equivalent to classical statistical mechanics, enabling theorists to apply a larger tool- 
kit to quantum field theory. For example, Wilson used a strong-coupling expansion to lowest 
order in 1/g^ to demonstrate confinement [4]. 

Among the techniques of statistical mechanics, the one that has become an industry is to 
integrate expressions of the form ([2]) on big computers with Monte Carlo methods. Although 
lattice gauge theory defines QCD mathematically and, thus, in principle provides an algorithm 
for computing anything, as with any numerical analysis, compromises are necessary in practice. 
To evaluate anything meaningful within a human lifetime, the integrals are defined at imaginary 
time, t = —ixi, turning Feynman's phase factor into the damped exponential of Eq. Q. A 
computer, obviously, has finite memory and processing power, so the spatial volume and time 
extent of the lattice are finite. 

From this expression, it is straightforward to derive some simple results for correlation 
functions. The two-point function 

(^(t)7rt(0)) = ^K0|7r|7r„)pexp(-m,„i), (3) 

n 

where vr is a composite field of definite quantum numbers (e.g., of the pion), and the sum ranges 
over all radial excitations. For time separation t large enough, a fit to an exponential yields 
the lowest-lying m^^ and |(0|7r|7ri)|. This is how we compute masses. For a transition with no 
hadrons in the final state, as in leptonic decays, simply replace 7r(i) with a current J: 

(7r(t)7rt(0)) = J^(0| JK„)(7r„|7rt|0) exp(-m^„t), (4) 

n 

in which the only new information is (0| J|7r„), yileding for large t the matrix element of the 
lowest-lying state. This is how we compute decay constants. For a transition with one hadron 
in the final state, one needs a three-point function: 

(7r(t)J(n)5t(o)) = ^(0|7r|7r„)(^,| J|B,„)(S^|5t|0) exp[-m^„(t - u) - ms^u] (5) 

mn 

in which the only new information is (7r„| J|i?m), yileding for large t and u the matrix element 
between the lowest-lying states. This is how we compute form factors and neutral hadron 
oscillation properties. 

Equations ©-([s]) are derived by inserting complete sets of eigenstates the QCD Hamiltonian. 
Thus, every successful fit of these formulae for hadronic correlators provides a posteriori incre- 
mental evidence that hadrons are indeed the eigenstates of QCD. 

In all cases of interest, the fermion action is of the form tpMip^ where the matrix M is 
a discretization of the Dirac operator (plus quark mass). Then the fermionic integration in 
Eq. pi) can be carried out by hand: 

(•) = ^ I VA [•'] det M exp (-fi-gauge) , (6) 

in which the fermionic integration replaces V'iV'i with [M~^]jj to yield •'. Importance sampling, 
which is crucial, becomes possible if det Mexp (— Sgauge) is positive. In most cases, a notable 
exception being the case of nonzero baryon chemical potential, this condition holds. 



The determinant det M is the mathematical representation of virtual quark-antiquark pairs, 
also called sea quarks. The matrix inverse Af~^ is the propagator of a valence quark moving 
through a stew of gluons A and sea quarks det M. Several quark propagators are sewn together 
to form hadronic correlation functions, from which masses and transition matrix elements can 
be computed. Computationally, detM is the biggest, and M~^ the second biggest, challenge 
in lattice QCD. The numerical algorithms become even more demanding as the quark mass is 
reduced, so in practice light quark masses are usually 2-5 times larger than the physical up and 
down quark masses. 

Because det M is so CPU- intensive, for many years lattice-QCD calculations were carried 
out in the "quenched approximation," in which detM is set to 1. Except for pilot studies of 
specialized new methods, the quenched approximation is now obsolete. It lives on in the jargon 
"unquenched lattice QCD," which means simply to do the right thing, namely compute det M. 
It also lives on in the jargon "partially quenched" QCD, which refers to unphysical set-ups in 
which Mvaience and Msea differ. (Often just the value of the quark mass differs.) This is useful 
because we believe we have a theory, a version of chiral perturbation theory, to incorporate the 
unphysical results into fits that, in the end, yield physical results [10] . 

How do the compromises of numerical lattice QCD impair the results discussed below? 
The imaginary time imposes no problem whatsoever for static quantities. The finite volume 
introduces errors that are exponentially suppressed and, hence, a minor source of uncertainty. 
Likewise for the finite time extent, except in thermodynamics (Sec. p|, where it becomes a tool. 
Finally the unphysical light quarks are exptrapolated away with self-consistent formulae from 
chiral perturbation theory [11]: a physical way to think about this step is that we remove the 
cloud of unphysically massive pions and replace them with the real thing (to some order in chiral 
perturbation theory). As in any industry, the techniques that work well for many things do not 
work for everything. 

By 2003 these techniques — including a realistic formulation of det M, had matured — making 
possible several postdictions [12]T. The next step was to make several predictions of hadronic 
properties that were not yet, but were soon to be, measured in experiments. These included 
form factors of D — t- Klu and D — ?■ iilv decays fr3|, the mass of the Be meson p^, and the 
decay constants of charmed mesons [T3]. For a description of these developments, see Ref. |16| . 

3. Hadron Spectrum 

We compute the masses of hadrons not only "because they are there," but also because it 
is interesting to see whether, and how, QCD generates mass. Particularly in extensions of 
the electroweak part of the Standard Model, there are many theoretical ideas for how mass is 
generated. Among these, however, we only know for sure that Nature makes use of the binding 
energy of gauge theories. Thus, QCD is a prototype for a mechanism that could reappear at 
shorter distance scales. 

It is instructive to start by thinking semiclassically about the potential energy V{r) between 
heavy quarks at separation r, or, equivalently, the force F{r) = —dV/dr. At short distances, V 
and F are Coulombic, with logarithmic corrections from asymptotic freedom. At large distances, 
the potential grows linearly, and, correspondingly, the force becomes a negative constant. An 
accurate picture of how this force arises is as follows. As r increases, a dipole field similar to 
that of electrodynamics forms. Gluons, being colored, attract each other, so the chromelectric 
field lines narrow, first into a sausage and eventually into a string. The QCD flux tube is full of 
energy, rising linearly with r. This energy, via m = E/c?, is the origin of hadron mass. 

This picture is confirmed in detail by lattice-QCD calculations of the potential, as shown in 
Fig. [TJ In addition to establishing the Coulomb+linear behavior described above (the points la- 
beled SjJ"), Fig.nishows an excitation of the chromoelectric flux between heavy quarks (the points 
labeled n„). This and further excitations elucidate how the QCD flux tube generates mass |17|y. 



In addition to this appealing picture, lattice 
QCD has been used to verify the mass 
spectrum of hadrons quantitatively, within 
a few percent. Figure [2] shows three sets 
of results, each with specific compelling 
features. Unless noted below, the error bars 
on the calculations encompass all systematic 
uncertainties. Figure [2^ shows the broadest 
attack on the spectrum [191 120]) including 
bb^ and cc states taken from Refs. [23l 123]. 
The baryon masses are less well determined 
than meson masses, partly because of larger 
statistical errors, and partly because the 
lattice formulation of light quarks in Refs. jl9| 
[2U] is sub-optimal for baryons. Figure ^p 
shows the spectrum with light in-simulation 
quark masses nearly as small as gl^^u + "i-d) 
[2T] . Such small quark masses are at the 
frontier, and the plot omits the error bar 
associated with the continuum limit. Finally, 
Fig. [2}; shows a complete calculation with good 
control of the baryons |22j . In particular, the 
nucleon mass, which provides almost all the 
mass in everyday objects [25], has now been 
verified within 3.5% to arise principally from 
chromodynamics: m = E/(?. 
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Figure 1. The heavy-quark potential V{r) 
vs. r, computed with lattice QCD (line and 
points labeled S+). From Ref. [E]. (The fine 
and points labeled II^j denote an excitation of 
the interquark chromoelectric field.) 
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Figure 2. Hadron spectrum computed with lattice QCD by (a) the MILC Collaboration [T91I20] . 
(b) the PACS-CS Collaboration [21], and (c) the BMW Collaboration [22]. 



4. Chiral Symmetry Breaking 

A striking feature of the hadron spectrum is that the pion has a smah mass, around 135-140 MeV, 
when most other hadrons have masses five or more times larger. For example, rrip = 770 MeV, 
nip = 938 MeV. To understand the origin of the difference, Nambu [6] applied lessons from 
superconductivity, noting (four years before quarks) that the pion's mass could be constrained 
to vanish by a spontaneously broken axial symmetry, with a small amount of explicit symmetry 
breaking allowing it to be nonzero. 

QCD explains the origin of this symmetry. If the up and down quarks can be neglected, the 
Lagrangian acquires an SUl(2) x SUr(2) symmetry, which provides a candidate axial symmetry. 
The consequences of spontaneous symmetry breaking were studied further by Goldstone [26], 
leading to a formula [27j , 

m'M^l^) = 0, (7) 

when applied to QCD with massless up and down quarks. The flavor-singlet expectation value 
{"tpi/j) is called the chiral condensate. If either factor on the left-hand side of Eq. ph is nonzero, 
the other must vanish. 

Twenty-five years ago, when the Lake Louise Winter Institute began, most physicists were 
confident that QCD was a good theory of the strong interactions, based on, for example, its 
explanation of Bjorken scaling in deep-inelastic scattering^. Because QCD was considered right, 
and because Nambu's picture of the pion was considered right, it was believed that QCD must 
generate a chiral condensate. There was, however, no direct calculation of {ipip) starting from 
the QCD Lagrangian, Eq. M. Now there is. Lattice QCD shows |28j 

(-0^) = [242 ± 4+11 MeV] ^ (MS scheme at 2 GeV), (8) 

where the first uncertainty is statistical and the second a combination of systematics. The 
quark masses have been adjusted to Nambu's idealization, rrin = nid — >• 0, rus physical. 
Spontaneous symmetry breaking also appears in two-fiavor QCD |29j . The chiral condensate 
has been established by direct computation to be very significantly nonzero. QCD breaks chiral 
symmetries spontaneously. 

5. SM Parameters 

The Standard Model has 19 free parameters, or 28 if nonzero neutrino masses and mixings are 
taken into account: 

• Gauge couplings: ctg, agED, "w = (A%/t')^/vr; 

• Lepton masses and mixing: m^^, m,,^^ ^^3, "^e, "i^, rrir] O12, 623, 6*13, (5pmns, "21, "31 j^ 

• Quark masses and mixing: muC*^, rrid, rrig, rric, m^, mt; Vijs, V^b> Kib> ^km; 

• Standard electroweak symmetry breaking: v = 246 GeV, A = (Mh/v)'^ /2. 

For ten or eleven of these parameters (bold), lattice QCD is either essential or important 
for determining their values of the natural world. Lattice field theory (without QCD) is also 
useful for shedding light on the Higgs self-coupling A [31j and the top-quark Yukawa coupling 
7/t = V^rrit/v [32], for which consistency of the field theory precludes arbitrary values. 

5.1. QCD Parameters 

Owing to confinement, there is no way to measure quark masses in a way comparable to, 
say, the electron mass. Instead, a Lagrangian definition as in Eq. (IT]), with suitable choice 
of renormalization scheme, must be determined from measurable properties of hadrons. For the 

^ For an explanation of neutrino mixing parameters, see Ref. |30j . 



Table 1. Quark masses from lattice QCD converted to the MS scheme and run to the scale 
indicated (first, second, third, and fifth columns). Also listed are charmed and bottom masses 
determined from e~^e~ production (fourth column). 





Ref. [20J 


Ref. [33J* 


Ref. IS] 


Ref. [35j 


Ref. [36]* 


fh^{2 GeV) MeV 


1.9±0.2 


2.01 ±0.14 


2.37±0.26 






?fid(2 GeV) [MeV] 


4.6 ±0.3 


4.79±0.16 


4.52 ±0.30 






ms(2 GeV) [MeV] 


88 ±5 


92.4 ±1.5 


97.7±6.0 






mc(3 GeV) [MeV] 








986 ± 13 


986 ± 10 


mb(10 GeV) [MeV] 








3610 ±16 


3617 ±25 



light quarks, the simplest hadronic property is simply the pseudoscalar mesons masses, with the 
physical value obtained when the pion and kaon masses agree with experiment. Three sets of 
results are shown in Table [I] The results in the second column |33j are derived from mass ratios 
underlying those in the first column |20j , as discussed below. The results in the third column are 
completely independent, in particular employing different methods for sea quarks and different 
approaches to electromagnetic effects. 

There are two noteworthy features of these results. First, the up and down masses are very 
small, about 4 and 9 times the tiny electron mass. Quark masses arise from interactions with 
the Higgs field, or its surrogate in other models of eletroweak symmetry. This sector is, thus, 
not the origin of much mass. Second, ?tt,u, though very small, is also very significantly far from 
zero. This is interesting, because were ttt-u = 0, then the additional symmetry of the Lagrangian 
would render 9 unphysical, obviating the strong CP problem. 

The heavy charmed, bottom, and top quark masses are large enough that they can be 
determined with perturbative QCD from features of high-energy scattering cross sections and 
energy distributions. For example, using perturbation theory to 0(ag) for the moments in s of 
the cross section for e~^e~ — )• cc, as a function of center-of-mass energy-squared s, one finds the 
result in the fourth column, fourth row of Table [1} The e'^e~ data can be replaced with moments 
of the charmonium correlation function, calculated with lattice QCD [371138]. Applying the same 
perturbative analysis yields the result in the fifth column, with astonishingly good agreement. 
The same methods can be applied to bottom quarks, also shown in Table [T} 

Returning to the light-quark masses, the new result of Ref. [33j is a precise value of the 
(scheme independent) ratio mc/rris = 11.85 ± 0.16. Combining this ratio with fhc |38jyand the 
ratios 2ms/(md ± m^) = 27.3 ± 0.3 and mu/md = 0.42 ± 0.04, both from Ref. |20], leads to the 
values in the second column of Table [TJ 

Lattice QCD also provides excellent ways to determine the gauge coupling Og = 5^/47r. In 
lattice gauge theory, the bare coupling Qq is an input. Alas, for most lattice gauge actions, 
perturbation theory in Qq converges poorly [39]"!', obstructing a perturbative conversion to the 
MS or other such schemes. Two other strategies are adopted to circumvent this obstacle. One is 
to compute a short-distance lattice quantity — a Wilson loop, a Creutz ratio, or the potential at 
separations of order a — and reexpress perturbation theory for the Monte Carlo results in a way 
that eliminates g^. The other is to compute a short-distance quantity with a continuum limit, 
and then apply continuum perturbation theory. The quarkonium correlator used for rric and mi, 
is an example: it also yields as{2mQ). Other examples include the Schrodinger functional |40j 
and the Adler function [41]. 

Results from several complementary lattice-QCD methods \38 \ l42fH3] are collected in Table^^ 
and compared to an average of determinations from high-energy scattering and decays [45j . One 
sees excellent consistency among results with different discretizations of the determinant for sea 
quarks. An important source of uncertainty is the truncation of perturbation theory, including 



Table 2. Values of as{Mz) from lattice QCD and an average of determinations from high- 
energy scattering and decays. A recent update to the values on the first two rows [38l HI] can be 
found in Ref. [30] . The central values and error bars from Refs. \i3\ I44j have been symmetrized 
to ease comparison. 

as{Mz) Observable Sea formulation Reference 

0.1183 ±0.0008 Wilson loops, Creutz ratios, etc. 2+1 asqtad staggered HPQCD \^^ 
0.1174 + 0.0012 charmonium correlator 2+1 asqtad staggered HPQCD |38j^^ 

0.1197 + 0.0013 Schrodinger functional 2+1 improved Wilson PACS-CS 

0.1185 + 0.0009 Adler function 2+1 overlap JLQCD f^ 

0.1186 + 0.0011 scattering, r decay, etc. 2+l(+l+l) Dirac (!) Bethke [45j 



strategies for matching to the MS scheme, and running to scale Mz- In the example of the 
lattice-scale loops, an independent analysis of the data from Ref. [l2] has been carried out, 
yielding OsiMz) = 0.1192 + 0.011 [l6]4^, to be compared with the first line of Table [2| 

As mentioned above, QCD is a union of the quark model of hadrons and the parton model of 
high-energy scattering. The agreement of the lattice-QCD results for Og, as well as for rric and 
?7ib, shows that hadrons and partons share the same QCD parameters, demonstrating QCD's 
breadth: the QCD of hadrons is the QCD of partons. 

5.2. Flavor Physics 

Like the quark masses, the quark mixing, or Cabibbo-Kobayashi-Maskawa (CKM), matrix |47| 
HH] arises from the electroweak interactions. In the Standard Model, the masses are proportional 
to the eigenvalues of the Yukawa-coupling matrices between the quarks and the Higgs doublet. 
The CKM matrix is the observable part of the transformations from the fields interacting with 
the weak gauge bosons to their mass eigenstates. Symmetries of the gauge interactions make 
many components of these transformations unobervable. For three generations, three mixing 
angles and one CP-violating phase remain to account for all the flavor- and CP-violation in 
nature. 

Lattice QCD calculations have played a key role in many aspects of flavor physics. A recent, 
comprehensive overview can be found in Ref. [l9j, so here we shall simply note some examples 
in which tension between the experiments and the Standard Model have appeared. These are 
tantalizing, because Standard CP violation seems insufficient to explain the baryon asymmetry 
of the universe. Tension appears in the global fit to the four CKM parameters [SO] and also in 
several specific flavor-changing processes. 

As anticipated, lattice-QCD calculations of neutral kaon mixing |5H [52] have improved such 
that the "standard" Standard-Model analysis had to re-incorporate certain effects of a few 
percent. (See Refs. |53l|5lj for details.) With the re-improved formula, the tension in the global 
fit strengthens |5ll [55] . 

Recent measurements of some purely leptonic decays B^ — )■ r^z^, Dg — ?• /U+z^, and Dg — )■ r^z^ 
are somewhat in excess of the Standard-Model prediction for the branching ratio. A crucial 
ingredient are the decay constants Jb and fos from lattice QCD [151 [56H58] . The nearly 2a 
excess of B'^ —)• r+z/ is usually interpreted as a possible signal of charged Higgs bosons [591160]. 
but then the non-Standard amplitude has to be around —110% of the Standard amplitude. 
The excesses of Dg — )• Z+z^ could be due to leptoquarks [61], with an few-percent amplitude 
constructively interfering. Note that the tension in this mode, which was once nearly 4(T, is now 
below 2(7, as improved calculations and, especially, new measurements have come out [62]. 



6. Thermodynamics 

Like any physical system, QCD has thermodynamic properties. During the early universe, the 
temperature was much hotter than it is now. In neutron stars, the baryon density is much 
higher than in normal nuclear matter. A sketch of the phase diagram is shown in Fig. [3^, based 
on lattice-QCD studies and models [63]. A surprising result is that the transition between the 
hadronic phase and the "quark-gluon plasma" is a smooth crossover, rather than a first- or 
second-order phase transition [65l [66]. This means that as the early universe cools, the hot 
matter becomes more and more like a gas of distinct hadrons. With a genuine phase transition, 
bubbles of the hadronic phase would form. At nonzero baryon density (chemical potential /u), 
it is thought that the transition becomes first order, but the matter is not yet settled ^^- For 
further discussion, see Refs \67\ 168], 

It may be worth clarifying what the quark-gluon plasma is. QCD thermodynamics is, as one 
would expect, based on the canonical ensemble, with thermal averages 



Tr 



(•) 



-H/T 



Tre-f^/^ 



(9) 



where T is the temperature. In quantum field theory formulated as in Eq. ([2]), the time extent 
A'^4 specifies a temperature T = (A^4a)~^. The trace Tr is over the Hilbert space of the QCD 
Hamiltonian H. The eigenstates — aka hadrons — do not change with T, but as T increases 
the propagation of a single source of color can change. First, thermal fluctuations encompass 
states with many overlapping hadrons, so color can propagate from one hadron to the next, 
as if deconfined. Second, the thermal average applies nearly equal weights to states of both 
parities, so chiral symmetry is restored — the thermal average of ipip vanishes even if the vacuum 
expectation value does not. With a smooth crossover, these changes need not emerge at the 
same T, but, in practice, it seems they do |65ll66j . The picture of thermal averages over hadronic 
eigenstates and the crossover nature of the transition may help us understand why hadron gas 
models of the transition are so successful. 

The nature of the QCD phase transition is influenced by the physical values of the light 
(up, down, and strange) quark masses, as sketched in Fig. |3J3. For vanishing quark masses, the 
transition would be flrst order. The ratio 2ms/(mu + ?7Zd) is well constrained by chiral symmetry 
(and substantiated by explicit calculation, as in Table [l]). But the masses are just large enough 
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Figure 3. QCD phase diagrams, (a) in the /i-T plane (image from Ref. [63]); (b) at ^ = 0, 
T = Tc, showing the order of the transition in the ?t1s-2('^u + 'n-d) plane (image from Ref. [64]). 



to push the QCD system into the region of crossover. If the Hght quark masses — cruciahy nis — 
were around half their physical size, the universe would cool through a first-order transition. 
What kind of fluke is this? 

7. Summary and Challenges 

During the twenty- five years of the Lake Louise Winter Institute, lattice gauge theory has 
developed several ideas about QCD, nurturing them from "QCD should work this way" to "QCD 
does work this way." Quantitative precision on a^ and heavy-quark masses reassures us that the 
QCD of hadrons is the QCD of partons. Accurate calculations of the hadron masses and the 
chiral condensate show how QCD generates mass and that QCD breaks chiral symmetry. The 
dependence of hadron masses on quark masses leads to the conclusion that the masses generated 
via interaction with the Higgs field is very small. Your mass is E/c? from QCD. Furthermore, 
the up quark's mass, though small, clearly does not vanish. The nonzero up, down, and strange 
masses make a qualitative difference to the early universe, because at physical quark mass (and 
low density) the QCD phase transition is actually a smooth crossover. 

On the quantitative front, QCD faces many challenges. Hints of non-Standard processes 
in flavor physics require ever more precise calculations. It is fairly certain that the Standard 
amount of CP violation is insufficient to explain the baryon asymmetry, so it is plausible that one 
of these hints will settle into real evidence. The advent of the LHC calls for other calculations 
that still lie beyond today's precision frontier of lattice QCD. For example, reliable moments of 
the gluon density inside the proton could help reduce uncertainties in LHC cross sections. If the 
LHC uncovers evidence for a dynamical mechanism breaking electroweak symmetry (similar in 
some, but not all, ways to QCD), lattice gauge theory will be necessary for non-QCD models 
(and, eventually, a new theory) [69l[70] . 

Nuclear physics is an arena where the need for computational lattice gauge theory is exploding 
[7T] . In many cases, the same basic methods [Eqs. ([3])-([5])] apply, but in others the technology has 
to be extended or invented [72j . Nuclear lattice QCD overlaps with astrophysics: better methods 
for nonzero baryon chemical potential would permit studies of the phases inside neutron stars, 
not to mention even denser phases shown in Fig. ^ ^73j; meanwhile, calculations of hyperon- 
nucleon interactions shed light on strangeness in neutron stars |74j . 
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